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Abstract 

High-resolution synchrotron x-ray diffraction measurements were performed on single crystalline 
and powder samples of B1M112O5. A linear temperature dependence of the unit cell volume was 
found between Tjy = 38 K and 100 K, suggesting that a low-energy lattice excitation may be 
responsible for the lattice expansion in this temperature range. Between T* ~ 65 K and Tjy, all 
lattice parameters showed incipient magnetoelastic effects, due to short-range spin correlations. An 
anisotropic strain along the a-direction was also observed below T* . Below T\r, a relatively large 
contraction of the a-parameter following the square of the average sublattice magnetization of Mn 
was found, indicating that a second-order spin hamiltonian accounts for the magnetic interactions 
along this direction. On the other hand, the more complex behaviors found for b and c suggest 
additional magnetic transitions below T/v and perhaps higher-order terms in the spin hamiltonian. 
Polycrystalline samples grown by distinct routes and with nearly homogeneous crystal structure 
above Tjv presented structural phase coexistence below T/v, indicating a close competition amongst 
distinct magnetostructural states in this compound. 

PACS numbers: 75.30.Kz,77.80.Bh,61.05.cp,61.50.Ks 
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I. INTRODUCTION 



Multiferroic materials with coexisting (anti)ferromagnetism and ferroelectricity have at- 
tracted renewed attention, due to the interesting physics involved as well as relevant potential 
applications in spintronics. A fairly strong coupling amongst ferroelectric and magnetic order 
parameters may occur as a result of exchange striction effects in magnetic structures lacking 
an inversion center (for a recent review, see Ref.-). Examples can be found within the class 
of frustrated antiferromagnets. The .RM^Os family is a particularly interesting case,- 1 ^ 
in which Mn 4+ C>6 octahedra and Mn 3+ Os pyramids are interconnected and no possible 
spin configuration can simultaneously satisfy all nearest-neighbor Mn-O-Mn superexchange 
interactions.-^ As a consequence, the magnetic structures actually found in this family 
frustrates some of the spin interactions. Below the magnetic ordering temperature, slight 
atomic displacements take place and strengthen (weaken) the satisfied (frutrated) interac- 
tions, breaking the inversion symmetry of the structure and leading to ferroelectricity.-^ 1 ^ 

While the above mechanism explains qualitatively the multiferroic behavior of .RM^Os, 
detailed experimental information on the atomic displacements associated with each of the 
commensurate or incommensurate spin structures of this family is still lacking. This is 
mostly due to the small magnitude of such displacements, presumably below ~ 0.01 A. 7 This 
limitation prevents a more quantitative test for the existing theories, most noticeably for 
ab-initio calculations. 7 On the other hand, the lattice parameters can be obtained directly by 
high-resolution x-ray diffraction experiments, and the thermal expansion coefficients may be 
also accurately obtained by macroscopic dilatometry measurements on single crystals. Such 
measurements may reveal the overall magnetoelastic coupling in the unit cell dimensions, 
and carry relevant information on the microscopic spin-lattice coupling mechanism leading 
to ferroelectricity. Dilatometry measurements of thermal expansion coefficients have been 
carried out for R = Ho, Dy, and Tb,- clearly revealing the lattice anomalies related to each 
of the magnetic transitions of these materials. 

It is well known that thermal expansion coefficients obtained by macroscopic dilatometry 
may have much higher resolution than the typical results of x-ray or neutron diffraction. On 
the other hand, diffraction is the only choice to investigate powder samples with anisotropic 
crystal structure, and may be also useful for single crystal studies if simultaneous lattice 
expansion and strain broadening measurements are required. This is generally the case 
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in investigations of ferroelectric materials, 9 which commonly have distinct behaviors for 
powder and single crystalline samples, also depending on details of sample growth. Thus, a 
technique combining the attributes of high resolution to detect minute lattice anomalies and 
microscopic sensitivity to probe powder samples and/or inhomogeneous phases is highly 
desirable to study ferroelectric materials, in particular the multiferroics. In fact, high- 
resolution synchrotron x-ray diffraction may be the technique of choice for some of these 
cases. 

BiMn 2 5 presents a magnetic structure with propagation vector r = (1/2, 0, 1/2) at low 
temperatures,— 1 ^ while other members show ground states with r = (k x , 0, k z ) with k x ~ 1/2 
and 0.25 < k z < 0,37.^&^2<^ Structurally, the major differences to the other compounds 
of the series arise from a largely distorted BiO§ cage, which has been ascribed to the electron 
lone pair in Bi 3+ .— No detailed study on the magnetoelastic properties of this compound has 
been performed, to our knowledge. To bridge this gap, we performed a synchrotron x-ray 
diffraction study on a single crystal of BiM^Os, as well as on two powder samples grown 
by distinct routes. It is shown that this technique may have enough resolution to reveal the 
relatively subtle lattice parameter anomalies related to the magnetic transitions and may 
provide quantitative details of the spin-lattice coupling in this compound. 

II. EXPERIMENTAL DETAILS 

The single crystal used in the present study were prepared by the flux method, as de- 
scribed elsewhere.- 1 ^ In addition, two polycrystalline samples of BMO were grown by entirely 
different routes. The sample named BMOl was grown by wet chemistry^ using the pro- 
cedure described in ref.— , while BM02 was grown by a solid state reaction: stoichiometric 
amounts of Bi2C>3 and Mn203 were thoroughly mixed in an agate mortar, and heated to 
800 °C for 48 h and to 900 °C for 72 h, with intermediate grindings at each 24 h. For 
powder diffraction experiments, a sieve was used to reject grains larger than ~ 5 /im, and 
the samples were deposited over flat Cu holders, apropriate for Bragg-Brentano geometry. 
Synchrotron x-ray diffraction experiments in the single crystal and powder samples were 
performed on the XPD beamline of the Laboratorio Nacional de Luz Smcrotron (LNLS),— 
using an incident beam with A = 1.5499 A for the single crystal study and A = 1.3773 A 
for the powder measurements, except when otherwise stated. The beam was focused onto 
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a spot ~ 0.8 x 2.0 mm 2 for the single crystal and ~ 0.8 x 5.0 mm 2 for powders at the 
sample position. A Ge(lll) analyzer crystal was placed in a goniometer attached to the 26 
arm, and a scintillation detector was used. The instrumental resolution for this setup was 
~ 0.01° full width at half maximum (FWHM) at 26 = 25°, 17 and a step width of 0.0025° 
was chosen for the powder measurements while the sample was rocked by 1° during each 
observation to minimize graininess effects. For single crystal measurements, a natural (001) 
surface with no treatment was chosen, and the lattice parameters were obtained from the 
Bragg positions of the (402), (045), and (006) reflections. Each reflection showed a clearly 
split three-peak structure in the axial (6 — 26) scans at all temperatures, indicating the 
presence of at least three domains with distinct sets of lattice parameters, possibly due to 
inhomogenously distributed Bi vacancies. All the domains showed identical temperature- 
dependence for the lattice parameters, and the results shown in this work (including the 
peak widths) are given for the strongest peak of each Bragg reflection. A closed cycle He 
cryostat was employed in our measurements, and the temperature was measured with an 
estimated accuracy better than ~ 1 K and stability of ~ 1 mK. All the measurements were 
performed below ~ 100 K, while the base temperature (11 K for powders and 17 K for 
the single crystal) was determined by the performance of the cryostat at the time of the 
experiments, dc-magnetic susceptibility measurements were performed on a commercial su- 
perconduction quantum interference device (SQuID) magnetometer, while the specific heat 
was measured on a commercial platform using the relaxation method. 

III. RESULTS AND ANALYSIS 

A. dc-magnetic susceptibility and Specific heat 

BiMn 2 05 has been described as a non-co linear commensurate antiferromagnet at low 
temperatures with propagation vector r = (5, 0, |), and the spins pointing nearly along the 
a-direction, as inferred from neutron powder diffraction (NPD) measurements by Munoz 
a/.— As a preliminary bulk characterization, we performed (ic-magnetic susceptibility (xdc) 
and specific heat (C p ) measurements on the single crystal (see Figs. [T|a) and (b)). Only 
one transition could be unequivocally distinguished in our C p data, within our resolution, 
while Xdc data show that the Mn spins are indeed oriented nearly along the a-direction in 
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the ordered phase, as inferred by the smaller susceptibility below T/v with the field along 
this direction. The Xdc data are consistent to those reported in Ref.—, and the conclusions 
drawn from Xdc and C p data are entirely consistent with the previous NPD results. 11 (xdc) 
measurements on the two powder samples also confirmed the AFM transition at T/v ~ 38—40 
K. 

B. X-ray diffraction 

Figures[2|a-c) show the temperature-dependence of the a, b, and c lattice parameters. The 
corresponding unit cell volume V is given in Fig. EJ^d). For temperatures between ~ T* = 65 
K and 100 K (our upper limit in this work), the evolution of all the lattice parameters follow a 
straight line within our resolution. The linear thermal expansion coefficients in this T-range 
are: a a = (l/a)5a/5T = 7.7(1) x lO" 6 K" 1 ; a b = 1.7(2) x lO" 6 K _1 ; a c = 3.27(4) x 10~ 6 
K -1 ; and the volumetric expansion coefficient is (3 = 1.27(3) x 10~ 5 K _1 . For T/v < T < T*, 
deviations from this behavior were observed. Interestingly, such deviations had a distinct 
sign for the a lattice parameter with respect to b and c, so that the unit cell volume followed 
the constant thermal expansion down to T/v ~ 38 K. At T/v, clear anomalies were observed for 
all lattice parameters and unit cell volume. While the a lattice parameter show a contraction 
below T/v, the b parameter show a peak-like feature close to T/v and c shows an expansion on 
cooling. The dominant contraction of a shown in Fig. Ufa) leads to a significant reduction 
of the unit cell volume below T/v (see Fig. 121(d)). 

In order to correlate the observed lattice parameter anomalies to the antiferromagnetic 
order parameter and obtain more detailed information on the magnetoelastic coupling in 
this material, the contributions from the non-magnetic thermal expansion coefficient were 
subtracted, leading to magnetoelastic contributions to a, b, and c, which we refer to as o^, 
and Cm- For this procedure, we assumed that the non-magnetic contributions to the 
lattice expansion have linear temperature-dependence also below T/v, at least down to 17 
K. These data were compared to the square of the average sublattice magnetization of the 
Mn 4+ (4/ site) and Mn 3+ (Ah site) ions (M 2 ), extracted from ref.—, also normalized at 17 
K (see Fig. |3]). It is interesting to note that the evolution of au closely follows M 2 below 
T/v- The peak-like feature in the bj^ parameter takes place at ~ 36 K, clearly below T/v, 
while Cm appears to show a feature at ~ 33 K, increasing linearly on further cooling. These 
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anomalies for &m and cm below Tjy are indicated by arrows in Fig. [21 

Figure H] shows the width of the (402) and (045) reflections, revealing a broadening of 
(402) on cooling below T* and a nearly constant width of (045). Such anisotropic strain 
broadening indicates a magnetically-driven fluctuation of the a lattice parameter throughout 
the sample, in contrast to b and c, within our resolution. 

Figure [5]^a) show the full powder diffraction pattern of BM02 at 100 K. The crystal 
structure was refined under the Pbam space group using the GSAS+EXPGUI suite.— ^ 
The structure reported in Ref.— at 300 K was used as the initial model for the refinement. 
The calculated profile after the refinement is also shown in Fig. [51 The experimental data 
with Q < 2.5 A -1 were excluded from this refinement to avoid the instrumental or extrinsic 
effects of peak asymmetry, self-absorption due to surface roughness and beam footprint larger 
than the sample size, which become noticeable at lower angles for the reflection geometry 
employed here. This procedure does not affect significantly the accuracy and precision of 
the refined parameters, since the density of Bragg peaks is larger in the higher Q-region. 
The overall fitting quality is satisfactory, and the refined structural parameters are given in 
Table [H No impurity phases were observed within our sensitivity. The equivalent pattern 
of BMOl (not shown) revealed weak unidentified impurity peaks (< 0.5 % of the strongest 
peak of the main phase. Overall, both samples were found to be homogeneous and of very 
good crystalline quality above T/y. No sign of anisotropic strain or symmetry lowering were 
observed for both samples at 100 K. 

Figures E](a) and[6](b) show a selected portion of the powder diffraction profiles for BMOl 
and BM02, respectively, including the (210) and (021) Bragg peaks, at several temperatures. 
It can be observed that, while the (021) peak shows nearly no T-dependence, the (210) 
reflection splits in two peaks at low T. A similar splitting was clearly identified in many 
other (hkl) reflections with a large h/(k + I) ratio. Attempts to index all Bragg peaks 
at 10 K within a single crystallographic phase with either a monoclinic or triclinic unit cell 
derived from the high-temperature orthorhombic cell were unsuccessfull. On the other hand, 
a model with two distinct phases with Pbam symmetry and slightly different sets of lattice 
parameters could match the observed Bragg peak splittings in the powder profile. These 
phases were labeled as PI and P2 in Figs. [6]^a) and[6](b). In order to avoid divergences in 
the fit, the b lattice parameter had to be constrained as equal in both phases. The refined 
lattice parameters for PI and P2 at 10 K are given in Table IIII1 It can be seen that PI has 



6 



smaller a and slightly larger c than P2. Unfortunately, the atomic parameters of PI and P2 
at 10 K could not be reliably extracted from the refinement, since the fit did not converge 
when all the relevant atomic parameters were simultaneously refined. This is mostly likely 
due to significant Bragg peak overlap of the two coexisting phases. 

To gain further insight into the nature of phases PI and P2 below T/v, the Bragg peak 
positions of the (210) reflections of both phases, obtained from a fit with using Lorentzian 
lineshapes (symbols), are given in Fig. Efa) for samples BMOl and BM02. For temperatures 
above T N , the positions of the single (210) peak corresponding to the unique structural phase 
are given. For the temperature interval near and below T/v in which the two peaks could not 
be reliably separated in the fit, the positions are not given in this figure. The "expected" 
positions of the (210) reflections, obtained from the single crystal data of Figs. [2](a-c), are 
given in Fig. [3(a) (solid line). An analysis of the (210) and (021) peaks widths obtained 
with a single Lorentzian fit for each reflection, is given in Fig. [T)(b) . While the (021) 
width shows only a weak temperature dependence, the (210) width increases steeply for 
temperatures below ~ T/v, for both BMOl and BM02 samples. This result shows that 
the onset of structural phase coexistence on powder samples takes place near the magnetic 
ordering temperature, T N ~ 38 — 40 K. 



IV. DISCUSSION 



We begin our discussion by considering the lattice parameter behavior above T/v (see Fig. 
[2]). In this temperature region, the thermal expansion is highly anisotropic, a a being much 
larger than on, and a c . Remarkably, the unit cell volume V shows a linear temperature depen- 
dence between ~ T/v and 100 K within our resolution, indicating a nearly constant thermal 
expansion coefficient (3. This is a highly unusual behavior, since f3 = b^T 3 + b^,T 5 + b7T 7 + ... is 
expected at sufficiently low temperatures,— and a constant (3 should occur only for T >> 9d, 
where Qd = 235 K is the Debye temperature for BiMn 2 05.— A possible explanation for this 
intriguing behavior is the hypothetical presence of a very low-frequency optical mode or some 
other dispersionless lattice excitation, which thermal population might dominate the lattice 
expansion below 100 K. This possible mode or excitation might be associated with coupled 
rigid rotations of the MnC>6 octahedra and MnOs pyramids, or most likely be a rattling mo- 
tion of Bi ions inside the BiOs cage. We should mention that an anomalous distribution of 
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Tb-0 distances was observed for the related compound TbM^Os, with relatively large ther- 
mal dependence even at low temperatures,— also suggesting low-energy lattice excitations in 
this family. Considering the linear volumetric thermal expansion above T/v, we expect that 
the energy of such mode should satisfy the relation E « A^T/v = 3.4 meV. Preliminary 
Raman^ and Infrared^ spectroscopy measurements did not reveal any clear optical phonon 
with energies between ~ 2 and ~ 5 meV, indicating that the suggested lattice excitation is 
either silent or have an energy below ~ 2 meV. We should mention that, since the linear 
temperature dependence of V occurs in the paramagnetic/paraelectric phase, it cannot be 

explained by the possible presence of electromagnons in the ferroelectric/antiferromagnetic 
phase.-2L28^ 

Even though the unit cell volume shows a linear temperature dependence down to T/v, 
the individual lattice parameters a, b, and c present deviations from a linear behavior below 
T* ~ 65 K. The sign of this deviation is opposite for a than for b and c (see Fig. [2](a-c)), 
compensating each other in the volumetric expansion between T/v and T*. Also, these rela- 
tively small deviations from the linear temperature dependence in the paramagnetic region 
have the same sign of the much larger magnetoelastic anomalies below T/v, i.e., a contraction 
for a and an expansion for b and c on cooling. This correspondence suggests that the devi- 
ations from linear temperature dependence between T/v and T* are not driven by phonons, 
but are rather related to the strong short-range spin correlations in this temperature region, 
consistent with a previous Raman scattering study in BiMn 2 05.— The magnetically-driven 
strain broadening of the (402) reflection below ~ T* (see Fig. HI) is an additional evidence 
that short-range spin correlations may influence the crystal lattice of BiMj^Os. 

The most noticeable lattice anomalies takes place below T/v, deserving a careful consider- 
ation. In a quadratic spin hamiltonian for transition-metal compounds, the magnetic energy 
and the atomic displacements due to exchange striction are proportional to the square of 
the magnetic order parameter (M 2 ), in a mean-field approximation, assuming that the an- 
gle between ordered spins remains constant below T/v. This proposition is valid for either 
Heisenberg, Ising, or Dzyaloshinskii-Moriya spin hamiltonians, or a combination of them. 
In order to accomodate the atomic movements directly related to exchange striction, the 
crystal lattice relaxes, possibly leading to complex structural changes even for atoms not 
directly related to the exchange mechanism. For small displacements, the magnitude of the 
elastic response of the lattice is directly proportional to the perturbing displacements due to 
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exchange striction, therefore the overal lattice anomalies arising from exchange striction are 
expected to follow M 2 for a quadratic spin hamiltonian. Figure [3]^b) shows that this simple 
prediction is confirmed for clm, but clearly fails for bu and cm in the whole studied temper- 
ature interval. In addition, these two parameters show interesting features at ~ 36 K and 
~ 33 K that might be associated to additional magnetic transitions below Tjy. We especu- 
late that these anomalies might be due to spin-flip and/or incommensurate-commensurate 
magnetic transitions, such as observed in other members of the family.- We should mention 
that no evidence of multiple magnetic transitions has been observed in our C p data, within 
our temperature resolution, while a two-peak structure in C p close to T/v was previously 
reported.— More detailed neutron diffraction experiments in the temperature interval close 
to T/v may be necessary to confirm or dismiss this hypothesis. 

The fact that only a,M scales with M 2 is interesting. This is also the direction in which 
the magnitude of the lattice anomaly is the largest. This can be rationalized on the basis of 
the magnetic structure of BiM^Os given in Ref.— . According to this, along the a-direction 
all the nearest-neighbor spin alignment between Mn 4+ 06 octahedra and Mn 3+ Os pyramids 
and between consecutive Mn 3+ Os pyramids in the sequence ...Mn 4+ -Mn 3+ -Mn 3+ -Mn 4+ ... 
are nearly AFM, with no clear manifestation of magnetic frustration except perhaps for a 
relatively small non-colinearity of the spin aligment. On the other hand, the spin alignment 
along b and c directions show clear signs of competing interactions. Along b, the coupling 
between Mn 4+ C>6 octahedra and Mn 3+ C>5 pyramids alternates between FM and AFM, break- 
ing the inversion symmetry of the structure and possibly causing ferroelectricity as in other 
members of the .RM^Os family. 1 Along c, the spin aligment between edge-shared Mn 4+ C>6 
octahedra also alternates between FM and AFM, the AFM and FM alignments correspond- 
ing to the shorter (=2.767 A) and longer (=2.988 A) Mn 4+ -Mn 4+ distances, respectively. 
Bi and Mn 3+ planes are intercalated between the AFM- and FM-coupled Mn 4+ planes, re- 
spectively. It has been argued that the Mn 4+ -0-Mn 3+ -0-Mn 4+ superexchange path and the 
longer Mn 4+ -Mn 4+ distances may help stabilizing the FM coupling among half of the Mn 4+ 
pairs. 11 We mention that, while the Mn spin structures in the ab plane are similar for all 
members of the family, the alignment along c is strongly dependent on the i?-ions, with the 
component of the spin propagation vector varying from k z = 1/4 for R =Er— to k z = 1/2 
for R =Bi.— >ii This is another clear indication of competing magnetic interactions along 
this direction. 
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The distinct features of the magnetic ordering along each direction obviously leads to 
an anisotropic magnetoelastic coupling, as manifested in the data of Fig. [2j First of all, 
the AFM alignment of nearest-neighbor spins along a leads to a relatively large contrac- 
tion of the lattice along this direction below T/v (see Fig. 12(a)), in order to enhance the 
AFM coupling of the ...Mn 4+ -Mn 3+ -Mn 3+ -Mn 4+ ... chains. On the other hand, some of the 
magnetic interactions are frustrated along b, and therefore the atomic distances related to 
'satisfied' pairs are expected to be reduced and others (between 'frustrated' pairs) should 
increase below T/v, leading to a compensation which is manifested as small values of bu in 
comparison to o« (see Fig. [3]). In the case of cm, a presumed increase of the separation of 
FM Mn 4+ pairs may be partially compensated by a corresponding approximation between 
the AFM Mn 4+ pairs or vice- versa, leading to smaller cm anomalies in comparison to au- 
In addition, the evolution of bu and cm may be also influenced by a lattice relaxation in re- 
sponse to the relatively large contraction of qm- Even considering such a complex situation, 
it is expected that magnetoelastic anomalies along all directions should follow the square 
of the sublattice magnetization for a spin hamiltonian composed only of quadratic terms 
on spin, as argued above. The fact that c and most notably b do not follow this behav- 
ior even for the temperature interval where no additional phase transitions are evidenced 
(T < 33 K) suggest that other terms, possibly of fourth order (biquadratic, three-spin and 
four-spin),— should be included in the spin hamiltonian to correctly account for the mag- 
netoelastic anomalies in BiM^Os. Although these terms are expected to be significantly 
smaller than the quadratic exchange terms, the magnetic frustration caused by the complex 
lattice geometry may cause a nearly complete cancellation of the quadratic terms, leading 
to a relative increase of importance of the quartic terms, and consequently a manifestation 
of such terms in the magnetoelastic anomalies in the directions where the spin frustration is 
most pronounced. This conclusion is likely extensible to the other members of the i?Mn 2 Os 
family and may be of relevance to a quantitative understanding of their complex magnetic 
phase diagram and multiferroic properties. 

The analysis of the powder samples provides relevant additional information. Contrary 
to the single crystal, a coexistence of two distinct phases with orthorhombic metrics (within 
our resolution) and slightly different sets of a and c lattice parameters was observed below 
~ T/v for both BMOl and BM02 powder samples. The lattice parameters of PI are closer to 
P2 in BMOl than in BM02, indicating that the magnitude of this effect in powder samples 
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is dependent on the sample details. The phase fractions of PI and P2 are also strongly 
sample-dependent (see Table HTT|) . Also, a comparison of the T-dependence of the position 
of a particular Bragg peak with temperature of both powder samples with that expected 
from single crystal data (see Fig. [7j) shows that the phase labeled PI (with smaller a lattice 
parameter) corresponds to the magnetostructural ground state of BiMn 2 Os, while P2 may 
be either a metastable state or a phase stabilized by slight chemical inhomogeneities (see 
below). It is well known that powder samples may favor the presence of metastable structural 
phases not generally observed on single crystals, leading to polymorphism. In the present 
case, the structures of PI and P2 converge above T/v, showing that the phase separation 
cannot be trivially explained by large chemical inhomogeneities between distinct grains. 
In fact, P2 presumably shows a different magnetic structure than PI, leading to distinct 
magnetoelastic coupling effects and therefore a slightly different set of lattice parameters 
at low temperatures. On the other hand, slight chemical inhomogeneities (such as in the 
Bi occupancy) not clearly evidenced in our powder diffraction profiles above T/v might still 
play a role in stabilizing a magnetically phase-separated state below T/v if there is a pre- 
existing close competition between distinct possible magnetic ground states. The fact that 
PI and P2 unit cells show differences in a and c lattice parameter with no observable change 
in b suggests that magnetic structures of these phases show distinct components of the 
propagation vector along a and/or c. Further studies are necessary to confirm or dismiss 
this hypothesis. 



V. CONCLUSIONS 



In summary, the magnetoelastic anomalies in single crystal and powdered BiMn 2 Os were 
investigated in detail by high-resolution synchrotron x-ray diffraction. It was found in the 
single crystal study that this compound shows a linear temperature dependence of the unit 
cell volume between T/v and 100 K, well below the Debye temperature, which has been 
associated to the possible presence of a low-energy lattice excitation that still remains to 
be directly observed. Below T/v, significant lattice parameter anomalies were found, most 
notably in the a-parameter, due to a magnetoelastic coupling. The anomaly of a was found 
to follow the square of the average sublattice magnetization, as expected for a quadratic 
spin hamiltonian. The anomalies of b and c follow a more complex behavior, signaling the 
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existence of additional transitions below T/v as in other members of the iflVh^Os family, 
and perhaps higher-order terms in the spin hamiltonian. It was argued that strong short- 
range spin correlations between T/v and T* = 65 K 21 give rise to an observable contribution 
to the linear thermal expansion coefficients. X-ray powder diffraction measurements taken 
on samples grown by distinct routes show phase coexistence of different magnetostructural 
states below Tjy. One of these states corresponds to that observed in the single crystal, 
while the other was attributed to a competing magnetic state which might be metastable or 
stabilized by subtle chemical inhomogeneities that might be present in the powder samples. 
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TABLE I: Refined lattice and atomic parameters of sample BM02 at 100 K. Errors in parentheses 
are statistical only, and represent one standard deviation. 



T = 100 K 


Pbam a = 


7.54116(1) A 


b = 8.52994(1) A c = 


5.75437(1) A 


Atom 


site 


X 


V 


z 


u iso (A 2 ) 


frac 


Bi 




0.15896(4) 


0.16556(4) 





0.00588(6) 


0.938(4) 


1V1111 


A f 


1/2 





n r ><\Qf\( r >\ 
u.zoyo^z ) 


u.uuzi J 


1 


Mn2 


Ah 


0.40755(15) 


0.35091(14) 


1/2 


.0029(2) 


1 


01 


Ae 








0.2876(10) 


0.0048(5) 


1 


02 


4<7 


0.1567(8) 


0.4453(6) 





0.0048(5) 


1 


03 


Ah 


0.1437(7) 


0.4243(6) 


1/2 


0.0048(5) 


1 


04 


8i 


0.3866(5) 


0.2018(4) 


0.2525(7) 


0.0048(5) 


1 


R p = 13.6 


% 




f^iup — 25.3 


% 




X 2 = 1-86 



TABLE II: Interatomic distances for sample BM02 at 100 K (A). Errors in parentheses are statis- 
tical only, and represent one standard deviation. 

Mn 4+ 6 Mn 3+ 5 

Mnl-02(x2) 1.961(4) Mn2-01(x2) 1.897(4) 

Mnl-03(x2) 1.872(4) Mn2-03(xl) 2.086(6) 

Mnl-04(x2) 1.923(4) Mn2-04(x2) 1.916(4) 

<Mnl-0> 1.919(2) <Mn2-0> 1.966(3) 

BiO s 

Bi-01(x2) 2.484(4) 

Bi-02(xl) 2.337(5) Bi-02(xl) 2.386(5) 

Bi-04(x2) 2.270(4) Bi-04(x2) 2.758(4) 
<Bi-0> 2.468(2) 
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TABLE III: Lattice parameters and unit cell volume for phases PI and P2 of powder samples 
BMOl and BM02 at 10 K. Errors in parentheses are statistical only, and represent one standard 
deviation. 



Sample 


Phase 


fraction (%) 


o(A) 


6(A) 


C (A) 


v(A 3 


BMOl 


PI 


90 


7.53039(3) 


8.52766(3) 


5.75426(3) 


369.517(3) 


BMOl 


P2 


10 


7.53730(7) 


8.52766(3) 


5.75377(10) 


369.827(6) 


BM02 


PI 


53 


7.53094(3) 


8.52955(3) 


5.75497(3) 


369.673(3) 


BM02 


P2 


47 


7.53888(4) 


8.52955(3) 


5.75405(3) 


370.004(3) 
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FIG. 1: (Color online) Temperature-dependence of (a) magnetic susceptibility taken with a mag- 
netic field of 2000 Oe along the a, b, and c directions; and (b) specific heat of the BiM^Os single 
crystal. 
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FIG. 2: (Color online) Temperature-dependence of the a, 6, and c lattice parameters, and unit cell 
volume V. Dashed lines fitting the data above 65 K are guides to the eyes. The vertical lines mark 
the magnetic ordering temperature TV and the onset of the strongly correlated paramagnetic state 
(see text and ref.— . 
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FIG. 3: (Color online) (a) Magnetoelastic component of a, b, and c lattice parameters, defined as 
the difference between the symbols and straight line in Fig. [2] and labeled as clm, ^m, and cm, 
respectively, (b) Scaling of au an d cm, normalized at 17 K (symbols), with the squared average 
sublattice magnetization of Mn + and M n 3+ ions (solid line) , extracted from ref.— . 
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FIG. 4: (Color online) Bragg peak widths obtained from radial (9 — 29) scans around the (402) 
and (045) reflections for the single crystal. The solid lines are guides to the eyes. 
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FIG. 5: (Color online) (a) Observed (symbols), calculated (line) and difference (botton line) x-ray 
powder diffraction profile of sample BM02 at 100 K, taken with A = 1.1271 A. The goodness-of-fit 
factors are, R p = 13.6 % and R p = 24.9 % (background subtracted), and x 2 = 1-85. The data 
were multiplied by five above Q = 5.3 A -1 for better visualization. The experimental data with 
Q < 2.5 A -1 were excluded from the refinement (see text). 
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FIG. 6: (Color online) (a) Small portion of the powder diffraction profile of BMOI and BM02 sam- 
ples covering the (201) and (021) reflections at selected temperatures. The profiles were vertically 
translated for better visualization. 
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FIG. 7: (Color online) (a) Symbols: temperature dependence of the positions of the (210) reflection 
for PI and P2 phases of BMOl and BM02 powder samples (see Fig. [6]). The solid line indicates the 
expected behavior taken from single crystal data shown in Fig. [2j (b) Temperature dependence of 
the width of the (210) and (021) reflections of BMOl and BM02, taken in the temperature region 
where the scattering at the (210) position could be fitted by a single lorentzian peak. 
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